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GLIMPSE-COl: the most massive intermediate-age stellar 
cluster in the Galaxy 



o 

(N 

> 

O 

O 



< 

o 

6 



> 

On 
(N 



13 



Ben Davies 1 ' 2 , Nate Bastian 3 ' 4 , Mark Gieles 4 , Anil C. Seth 5 , Sabine Mengel 6 , and 
Iraklis S. Konstantopoulos 7 

1 Rochester Institute of Technology, 54 Lomb Memorial Drive, Rochester, NY 14-623, USA 

2 School of Physics & Astronomy, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK 
^School of Physics, University of Exeter, Stocker Road, Exeter EX 4 4QL, UK 

4 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK 
5 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA 
6 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 

7 Eberly College of Science, The Pennsylvania State University, 525 Davey Lab, University Park, PA 16802, USA 



ABSTRACT 

The stellar cluster GLIMPSE-COl is a dense stellar system located in the Galactic 
Plane. Though often referred to in the literature as an old globular cluster traversing 
the Galactic disk, previous observations do not rule out that it is an intermediate 
age (less than a few Gyr) disk-borne cluster. Here, we present high-resolution near- 
infrared spectroscopy of over 50 stars in the cluster. We find an average radial velocity 
is consistent with being part of the disk, and determine the cluster's dynamical mass 
to be (8 ± 3)xl0 4 M©. Analysis of the cluster's M/L ratio, the location of the Red 
Clump, and an extremely high stellar density, all suggest an age of 400-800Myr for 
GLIMPSE-COl, much lower than for a typical globular cluster. This evidence therefore 
leads us to conclude that GLIMPSE-COl is part of the disk population, and is the 
most massive Galactic intermediate-age cluster discovered to date. 
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1 INTRODUCTION 

With the advent of large near/mid-infrared (IR) detec- 
tors and wide-area surveys, many massive stellar clusters 
(> 10 4 Mq) have recently been discovered in the plane of 
the Galaxy. So far, due to the selection techniques, only 
extr emely young clus t ers (< 25 Myr) have been identife d 
(e.g. IClark et al.ll2005l : iFiger et al.ll2009 : iDavies et al.1l2007h 
Howe ver, based on size-of-sample effects ( Gieles fc Bastianl 
l2008f) and what we see in similar face on spiral galaxies 
( Larsen fc Richtlerll 19991 ) we expect a large number of inter- 
mediate age (100-1000 Myr) clusters with masses between 
(1O 5 -1O 6 M0). Due to their stellar population properties, 
these clusters are expected to be very difficult to detect due 
to confusion with background stars. 

A candidate for such an intermediate age cluster may 
be GLIMPSE-COl (hereafter GC01), an object discovered in 
the Spitzer GLIMPSE m id-IR survey of of the Galactic disk 
l)Kobulnickv et aD 120051 . hereafter K05) . Despite the clus- 
ter being located 0.1 degrees from the Galactic plane, these 
authors suggest that G01 is an old globular cluster which 
happens to be passing through the disk of the galaxy. How- 
ever, they note that the colour-magnitude diagram (CMD) 



is equally well fit by a intermediate age stellar population 
(< few Gyr; see Fig. 7 in K05), while it may be argued that 
the cluster's luminosity function provided inconclusive re- 
sults (their Fig. 8). Further studies of the cluster have not 
pro vided convinc i ng evi dence of the cluster's age or metallic- 
ity llvanov et al.l (|2005f ) provided higher resolution near-IR 
photometry, detecting what appears to be the Red-Clump in 
the colour-magnitude diagram (see their Fig. 4) , and argued 
from the slope of the Red Giant Branch (RGB) that the 
cluster was metal poor. However, as stated by Ivanov et al., 
their method makes the implicit assumption that the cluster 
was >10Gyr old, and is not calibrated for intermediate-age 
clusters. 

In this paper we present high resolution near-IR spec- 
troscopy of several stars in the cluster. From these data we 
are able to determine radial velocities for each star observed, 
measure the cluster's velocity dispersion, and ultimately the 
cluster's dynamical mass. Additionally, we obtain a very re- 
liable measurement of the cluster's average radial velocity, 
which we can compare to the Galactic velocity field to estab- 
lish whether or not the cluster is co-moving with the disk. 
We combine this new information with that already present 
in the literature to attempt to shed new light on the nature 
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Figure 1. Examples of the data obtained, while in red we show 
the spectra prior to the fringing correction. In order from top to 
bottom, wc show a bright object (K=8.31) with strong CO; a 
bright object (A"=9.68) with weak CO; a faint object (K =10.87) 
with strong CO; and finally the spectrum of the faintest star we 
obtained (fsT=12.13). 

of GC01 - is it an old globular cluster passing through the 
disk, or is it an intermediate age disk-borne cluster, the first 
of it's kind detected in our Galaxy? 

We begin in Sect. [2] with a description of the observa- 
tions, data reduction procedures and analysis techniques. In 
Sect. [3] we describe our results, and discuss the issues of the 
cluster's physical properties such as distance and mass. The 
nature of GC01 is discussed in Sect.[4l while we conclude in 
Sect. 



2 OBSERVATIONS & DATA REDUCTION 

Our data were taken as part of the ESO observing pro- 
gramme 383.D-0025 (PI N. Bastian). We obtained near- 
infrared spect ra of numerous stars in the field of GCOl 
using ISAAC (|Moorwood et alJll99sh . mounted on UT1 of 
the ESO-VLT. We used the instrument in medium resolu- 
tion mode with the 0.3" long-slit and central wavelength of 
2.3/im. This set-up achieves a spectral resolution of i?=8900 
in the range ~2.25-2.35/xm, allowing us to observe the CO 
bandhead at 2.293/im and neighbouring continuum. 

We employed the observing strategy of stepping the slit 
across the cluster to positions aligned with bright stars. This 
technique also obtains data on fainter stars which fall into 
the slit serendipitously at each position. The observations 
were split across 2 observing blocks (OBs): one which in- 
cluded the seven slit positions covering the brightest stars, 
with integration times of 60s (NDIT=1); and the remaining 
30 slit positions which had 90s (NDIT=1) integrations. At 
the end of each OB arc lamp observations were taken for 
wavelength calibration purposes. We observed the B5v star 
Hip094378 to measure the telluric absorption immediately 
after the science observations. Flat-field and dark frames 



were taken at the end of the night, while at the beginning 
of the night a bright standard star was observed at several 
stepped positions along the slit to characterize any spatial 
distortion on the detector. 

Our data reduction began with subtraction of the dark 
current and dividing though by a normalized flat field. The 
spatial distortion of the 2-D frames was determined from 
the stepped standard star observations, while the distortion 
in the dispersion direction was measured from the arc-lamp 
observations of each observing block. The distortion in the 2- 
D frames was then corrected by resampling each frame onto 
an orthogonal grid, resulting in an absolute wavelength cal- 
ibration accurate to ±2 km s _1 (as measured from the resid- 
uals between the observed and literature values of the arc 
lines). The sky emission in each 2-D frame was character- 
ized by measuring the background in regions containing no 
star traces and interpolating across the full length of the slit. 
Stellar spectra were extracted by summing across the rows 
at each channel, being careful to avoid contamination by 
neighbouring stars in crowded regions. Telluric absorption 
was removed by dividing through by the standard star spec- 
trum, after first correcting for the star's intrinsic continuum 
slope by dividing through by a black-body curve. 

In the reduced spectra it was apparent that there was 
a degree of fringing of up to 10% of the continuum, towards 
the red end of each spectrum, with a period of ~0.04/im. 
To correct for this we identified a spectrum in each observ- 
ing block that had no detectable spectral features (i.e. no 
CO bandhead) and had high signal-to-noise. These stars are 
likely to be early ~G-type stars. This spectrum was then 
smoothed with a narrow filter to remove noise and intrinsic 
features, and was then used as a measure of the fringing pat- 
tern. The fringing pattern of each target star is a function 
of the star's position on the slit, so the 'fringe' spectrum 
was first cross-correlated with each science spectrum before 
dividing through. In Fig. Q] we plot four example spectra, 
showing the quality of our data. We show bright and faint 
stars, as well as stars with both weak and strong CO ab- 
sorption. The spectra prior to the fringing correction are 
also shown to illustrate that the magnitude of this artifact 
is only minor. 

To measure the radial velocity of each star observed, we 
cross-correlated each sp ectrum with that of Arcturus (from 
I Wallace fc Hinklelll996l ), after first degrading the Arcturus 
spectrum to match the spectral resolution of our obser- 
vations. This cross-correlation procedure provides relative 
velocities accurate to a few one-hundredths of a pixel, or 
<1 km s _1 . At this level we are likely dominated by system- 
atic errors, so for a conservative estimate of the uncertainty 
on each velocity measurement we use the absolute error in 
the wavelength calibration solution, ±2 km s _1 . Compar- 
ing the radial velocities of the stars obtained from each of 
the two observing blocks, we find that the two averages are 
within ±1.5 km s _1 of one another, indicative of the absolute 
error on the wavelength calibration. 



3 RESULTS 

In Table Q] we list the observed properties for all stars ob- 
served for which we could obtain corresponding near-IR pho- 
tometry. Astrometry for each star was determined from the 
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Table 1. Observational data of the stars observed. Astrometry is taken from slit positions, and is accurate to ±1". EW(CO) was measured 
between 2.293^tm and 2.310/im, uncertainties are ±5%, estimated from repeated measurements. E rrors on veloc i ties a re conservatively 
±2 km s _1 , dominated by the uncertainty in the wavelength calibration. Photometry is taken from llvanov et all l|2005h . 
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Figure 2. Left: A A*-band wide-field image of the cluster taken from UKIDSS, with the stars for which we obtained spectroscopy 
identified by circles. Right: An illustration of the cluster's velocity dispersion. Symbol sizes denote the magnitude of the stellar velocities, 
once the average cluster velocity has been subtracted. Objects with relative velocities greater than 30 km s — 1 were deemed to be unrelated 
foreground objects, and are identified by crosses. 
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Figure 3. Histogram of the observed velocities. The mean veloc- 
ity and its standard deviation (clipped at 92±30km s — 1 to pre- 
clude foreground objects) is illustrated by the blue long-dashed 
line. The red dotted and dashed lines show the velocity distribu- 
tions of the two separate observing blocks in which the data was 
taken. 



telescope pointing and the star's position on the slit, with 
UKIDSS images used to fine-tune the absolute calibration 
of the RA and DEC in the file headers. Photometry was ob- 
ta ined by cross-correl ating with the photometric catalogue 
of llvanov etaD (|2005h . For six stars in our sample no cor- 
responding photometry could be found. These data were 
discarded from the sample, since the signal-to-noise ratio 
(SNR) was much lower than the rest, which typically have 
SNR>50. We checked for systematic errors in our data by 
checking for correlations between radial velocity and bright- 
ness, colour, slit position, coordinates and CO absorption 
strength, finding no apparent trends. 



3.1 Radial velocities 

The stars observed are indicated in Fig. [2] while in the right- 
hand panel of the figure we illustrate their spatial velocity 
distribution. A histogram of the observed velocities (in the 
Local Standard of Rest frame) is shown in Fig. [3] The major- 
ity of the velocities are within 90±30km s _1 , with a small 
number of measurements of much lower velocities. These 
objects with low line-of-sight velocities are likely to be fore- 
ground objects which are not physically related to the clus- 
ter. To find the average velocity of the observed stars we 
took the mean of all measurements iteratively clipped at 3<r 
to discard objects with outlying velocities. We find a mean 
velocity of 91.6km s _1 , with clipping limits ±30 km s _1 of 
this value, and 53 stars having velocities falling within this 
range. The formal error on the mean velocity (= a-o/s/n) is 
compounded by the absolute error on the wavelength cal- 
ibration, therefore our formal measurement of the average 
cluster velocity is 91.6 ± 2.4 km s _1 . Altering the clipping 
limits by 10 km s" 1 in either direction does not affect the 
average velocity outside this error margin. 

To find the velocity dispersion, we calculate the r.m.s. 
standard deviation on the sigma-clipped mean, and subtract 
the absolute error in quadrature. Before making this calcu- 
lation however, we first assume that the stars observed in 
each OB are sampled from the same velocity distribution, 
shifting each population by ±1.8 km s _1 so that they have 
the same mean. Following this procedure, and clipping at 
92±30 km s _1 , we find a v — 9.6 km s _1 . If we do not shift the 
velocities of the two OBs, we find a v = 9.9kms _1 . If we nar- 
row the clipping limits by 5 km s _1 wc find a v = 9.0 km s _1 . 
Hence, our measurement is stable to within ±0.6 km s _1 . 

To estimate the formal empirical error on a v , we ran 
monte-carlo experiments in which we generated a ran- 
dom set of velocities with a gaussian distribution and 
rj=9.6km s _1 , and folded in the experimental uncertainty 
of ±2 km s _1 . We then determined what would be the ob- 
served velocity dispersion once the instrumental resolution 
was removed. We ran these experiments for several sizes of 
samples, from 10 to 100 measurements, and calculated the 
empirical error for each. We found that for a sample size of 
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consistently 1.1km s . Our formal measurement of GCOl's 
velocity dispersion is therefore a v = 9.6 ± 1.1 km s~ . 



3.2 Cluster rotation 

In addition to the velocity dispersion, we also investigated 
the cluster for evidence of rotation. The relaxation time for 
a cluster of GCOl's mass and size (see later) is ~ 10 8 yrs, so 
it is unlikely that any rotation would be reta ined if it had 
an ag e typical of globular clusters, ~10Gyr l|Fall fc Frenkl 
1 19851 ). Instead, the presence of rotation would suggest a 
much younger age. 

We divided the cluster into two with a line going 
through the cluster centre which had a position-angle 8. We 
then considered the mean velocities of the stars on oppo- 
site sides of this line as a function of 9. We found that the 
maximum velocity difference was ~2.5km s _1 at an angle 
of 6 = 130 ± 40° east of north (indicated on the right panel 
of Fig.0. This is consistent with being perpendicular to the 
cluster's apparent elongation (~60°, see K0fJ3), suggesting 
an oblate-spheroidal structure to the cluster. 

To test the significance of this result we performed a 
number of trials with randomized the velocities, keeping the 
stellar positions fixed and repeating the same test for ro- 
tation. We found that our result of a velocity gradient is 
significant at the 1.8<r level. Hence, any evidence we see for 
rotation is marginal. We also note that any rotational veloc- 
ity contribution appears to be much smaller than the overall 
velocity dispersion. 



3.3 Distance 

K05 use the 13 CO integrated column density along the line- 
of-sight to GC01 to constrain the distance to the cluster. 
They argue that the extinction to the cluster is consistent 
with being behind a 13 CO cloud at «LSR=46km s _1 , but in 
front of clouds at 81 km s _1 and 100 km s _1 . As the latter ve- 
locity corresponds to the tangent point at 7.3kpc, while the 
fLSR=46km s -1 cloud has a kinematic distance of 3.1kpc, 
K05 conclude that the distance to GC01 must be 3.1-7.3kpc. 

If, for the moment, we assume that the cluster is part 
of the disk rather than a globular passing through the 
disk, then the measured velocity «LSR=90±4km s _1 yields 
kinematic distances of d nC ar = 5.0 ± 0.9k pc and rif ar = 
7.9 ± 0.4k pc (using th e rotat ion curve of l|Brand fc Blitj 
Il993h : see iDavies et all (|2007l) for more details). If we dis- 
card the far-side velocity on the basis that the cluster is 
unlikely to be beyond the tangent point (K05), then we 
find a consistent picture whereby the cluster is located at 
5.0±0.9kpc, behind the molecular cloud at 3.1kpc which is 
responsible for the majority of the line-of-sight extinction. 
Thus, GC01 radial velocity is consistent with the object co- 
moving with the disk. 

In another study of the cluster's distance, llvanov et al.l 
|2005h detected what appears to be the Red Clump in a 
near-IR photometric study of GC01, located at K ~13 



1 The Spitzer/GLIMPSE images of GC01 indicate that the clus- 
ter suffers from differential reddening, which may give the ap- 
pearance of an elongated morphology 
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Figure 4. Hess diagram for two synthesized c lusters, constructed 
using the evolutionary models of ICioni et all ll2006al rbT). The plot 
shows a metal-rich lGyr-old cluster (black solid contours) and 
a metal-poor, lOGyr-old cluster (blue dashed contours). In the 
regime of current observations (K <14), the clusters are indis- 
tinguishable. Both have stellar overdensities at X=13, and have 
RGB slopes that may appear similar if the upper RGB is poorly 
populated. 
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Figure 5. The distance to GC01, inferred from the IR bright- 
ness and colour of the Red Clump, depending on the c luster's 
age, using the evolutionary models of lCioni et a l. (2006a b). The 
upper and lower limits to the cluster's distance, argued by K05 
on the basis of the 13 CO column density, are marked by the green 
dashed lines. 
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Figure 6. Left: Surface-brightness profile of the cluster from observations at different telescopes and through different filters, all of which 
give very similar values for R c g (indicated by the dashed line). The shaded area indicates the magnitude of the uncertainty on i? e g. 
Right: UKIDSS _R"-band image of the cluster, logarithmically scaled between the background level and IOOOct above the background, with 
the half-light radius illustrated by the dashed circle. 



and J — K ~3. These authors used this feature to esti- 
mate the distance to the cluster. However, they only used 
a single calibration for the absolute magnitude of the Red 
Clump, when in fact the location of this feature depend s 
upon cluster age and metallicity (Pie trzvnski et all [2003). 
Here, we have re-evaluated the distance to GC01 using 
the Red Clump me thod. Taking the evolutionary models of 
ICioni et al.l J2006allbl . and references therein) and a Kroupa 
IMF ( Kroupall2002r i. we generated synthetic populations of 
stars for a range of ages and in two distinct metallicity 
regimes: Solar (i.e. metal-rich) and l/20th Solar (i.e. metal- 
poor). For each synthetic cluster we created a K vs J — K 
colour- magnitude diagram (CMD), and adjusted the clus- 
ter's extinction and distance until a good match to the sim- 
ulated and observed Red Clump was found. We show two 
such simulations in Fig. U where it can be seen that clusters 
with very different ages and metallicities can produce stellar 
overdensities at the same location in the CMD. 

The results of this experiment are shown in Fig. [S] The 
age and metallicity dependence on the Red Clump distance 
calibration is rather weak for clusters older than ~500Myr, 
while for clusters younger than this the relation is very sen- 
sitive to age. We see that, within the distance upper and 
lower limits argued by K05, values of log(age/Myr) = 8.5 
- 10 (300Myr - lOGyr) are permitted at either metallicity. 
The implied extinction is also dependent on the cluster age, 
ranging from Ak = 1.52 — 1.61 in the metal-rich models and 
Ak = 1-63 — 1.73 in the metal-poor models. 

To summarize, we have insufficient evidence to point 
towards a definitive cluster distance within the upper and 
lower limits of 3.1-7.3kpc. If the cluster is comoving with 
the Galactic disk, as suggested by its radial velocity, then 
the cluster is unlikely to be older than ~lGyr, and so we 
have a consistent picture from the Red Clump and kine- 
matic distances in which the cluster is located at «5kpc. 
However, if the cluster is an old globular, with a low metal- 
licity and an age J>10Gyr, then the location of the Red 



Clump suggests that the cluster distance is close to the lower 
li mit ( 3-3. 5kpc) , cons istent with the distance determination 
of llvanov etafl (120051 ). 3.7±0.8kpc. In this case the cluster 
would likely be traversing the Galactic disk, and the simi- 
larity between the cluster's kinematic and physical distances 
would be purely coincidental. Comparing the observed ra- 
dial velocity with the distribution of old globular clusters, 
we find that there is a 10% chance that it has a velocity 
that happens to match the expected galactic rotation, cal- 
cu lated using t he velocity distribution of globular clusters 
m lHarrisI l| 19961 ). Hence, though the radial velocity suggests 
that GC01 is part of the disk population, this evidence alone 
cannot provide definitive evidence one way or the other to 
the cluster's nature. 



3.4 Cluster size 

In order to estimate the effective radius R e g of the clus- 
ter we utilise imagin g from the 2MASS (J, H, and Ks) 
(|Skrutskie et all 120061 ) and UKIDSS (we only used the J- 
band images as the H and K-band images suffered severe 
saturation effects) near-IR surveys. In both cases we esti- 
mated the position of the cluster center using isophotes at 
distances larger than 10" from the determined center. We 
then summed the flux in concentric circles, centered on this 
point, and subtracted the (area normalised) flux in a back- 
ground annulus > 40" from the cluster center. This cumula- 
tive flux distribution was then normalized to the maximum 
value attained, and the effective radius was defined as the 
radius containing half the light of the cluster. Our measure- 
ment of the cluster's flux profile is shown in the left panel 
of Fig. El 

In order to test the stability and estimate the error of 
our effective radius determination, we carried out the same 
analysis adopting different cluster centres in a box of width 
and height of 8" centred on the original position. Addition- 
ally, we varied the background annulus from 40" to 90". 
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The mean of these tests was _R e ff = 14.5". This is shown as a 
dashed circle in the right panel of Fig. [6] The error was then 
estimated as the standard deviation of all of these trials, and 
found to be 3". Overall, the agreement between the derived 
effective radius using the 2MASS and UKIDSS images was 
very good. 

This radius is significantly smaller than that reported 
by K05 (36"). As an extra check on our radius measure- 
ment, we repeated our analysis on the GLIMPSE images of 
the cluster, finding measurements consistent with that from 
near-IR images (14.1±0.8" and 13.8±0.8" for the 3.6A«n and 
4.5/im images respectively). It is not clear where the differ- 
ence between our measurement and that of K05 stems from, 
but one possibility lies in the estimate of the (high) back- 
ground, due to the location of the cluster in the Galactic 
plane. 

Assuming a distance to the cluster of 3.1, 5.0 or 7.3 kpc 
(see Section 3.1), results in an effective radius of 0.22, 0.35 
and 0.51 pc, respectively. Hence, G01 is a very dense cluster 
and app ears to be mu ch more compact than typical globular 
clusters i|Harrisll 19961 ) and yo ung/intermediate aged m assive 
clusters in external galaxies (|Scheepmaker et al.ll2007T ). 

3.5 Cluster mass 

We now calculate the mass of GC01 under the assumption 
that the cluster is in virial equilibrium. It has recently been 
argued that for very young clusters this assumption may 
not be valid, as the removal of intra-cluster gas by super- 
novae leaves the velocity dispersion of the cluster super- 
virial. However, this effect is tho ught to be important only 
for clusters with ages <2 0Myr l|Bastian fc Goodwin! 120061 ; 
iGoodwin fc Bas tian 2006). Binaries can also influence the 
velocity dispersion of a stellar cluster, though this effect be- 
comes less important for clusters older than ~100Myr, es- 
peciall y if the cluster is massive and compact (IGieles et al.l 
l2010bl ). We neglect any contribution from binary motions in 
our analysis of the velocity dispersion. 

To determine the dynamical mass Md yn we use the 
equation, 

Afdyn = T) G (1) 

where 7? c fr is the cluster radius, G is the gravitational con- 
stant, and r\ depends on the cluster density distribution but 
is typically taken to be «10. Using the value of R c \ = 14" 
measured earlier, we find a cluster mass of (8 ±3) x 10 4 M o , 
where the error is dominated by the error in distance. While 
this is a typical mass for a globular cluster, young clusters 
in this mass range have also been observed in the Galactic 
plane (e.g. Wd 1, RSGC2 lBrandner et al Jl2008l ; lDavies et all 
120071 ). Indeed, for a disk-borne cluster to survive until an age 
of ~lGyr without bei ng dissolved it wou ld likely need to be 
at least this massive l|Lamers et al.|[2005l 'l. 

GC01 appears to be an extremely dense and compact 
cluster, with a density of \og(ph / Mq pc~ 3 ) = 5.0±0.4 (where 
ph = 3M c i/(87ri?J), and the half-mass radius Ru = 4/3-Reflf). 
Here, the uncertainty is again dominated by that on the 
cluster's distance. This is de nsity exceeds that any globu- 
lar cluster in the catalogue o f lHarris! l|l996h : the most dense 
is NGC 6540 with \og{p h ) = 4.9, while all globular clus- 
ters with masses similar to GC01 have densities an order of 



magnitude lower 0. The density is however co mparable to 
very young mass ive clusters, such as the Arches (|Kim et al.l 
l2006h . Recently, IGieles et al.l (|2010aT ) have performed ana- 
lytic and N-body calculations to look at the evolution of 
stellar clusters, in particular their size and density. These au- 
thors find that clusters cannot exist for indefinite periods in 
extremely dense state, but rather expand on the timescales 
of a few relaxation times. If GC01 is an old globular, we 
would have expected it to have expanded during the past 
10-12 Gyr, and hence the current dense state of GC01 ar- 
gues for a much younger (<Cl Gyr) age. 



4 DISCUSSION: THE NATURE OF 
GLIMPSE-COl 

We now address the topic of GCOl's nature - specifically, is 
it an old, metal-poor globular cluster passing through the 
plane of the Galactic disk, or is it a young, metal-rich disk- 
borne cluster? 

The discussion in K05 provides convincing arguments 
against the cluster being very young, i.e. <50Myr, due to the 
lack of bright Red Supergiants, OB stars and intra-cluster 
material. However, as we have already discussed in Sect. (TJ 
their analysis is unable to distinguish conclusively between 
the cases of a very old metal-poor (i.e. globular) cluster and 
that of a metal-rich intermediate age (~lGyr) cluster. 

We also mention t he detect i on of extended X-ray emis- 
sion around GC01 by iMiraball (|2010l ) , spatially coin cident 
with the mid-IR emission seen by K05. IMiraball (|2010l ) spec- 
ulate that this emission may be arising in a bow-shock struc- 
ture, created as the cluster plunges through the disk, al- 
though they note that there are other explanations for the 
emission which cannot be ruled out. 

In this paper we have shown from GCOl's velocity dis- 
persion that the cluster's mass is (8 ±3) x10 4 Mq. Th is mass 
is typical for globular clusters (see e.g. lHarrislll999l '). while 
it is certainly not unheard of for young Galactic clusters ei- 
ther. Hence, the cluster mass does not provide a persuasive 
argument for either side of the debate. 

Similarly, the cluster's radial velocity does not rule out 
either argument. A radial velocity inconsistent with the local 
Galactic rotational velocity would have been a strong indi- 
cation that the cluster was not co-moving with the disk, and 
so was likely an old globular passing through. However, the 
velocity we measure is consistent with the Galactic rotation 
curve at a distance of ~5kpc, which is perfectly within the 
upper and lower distance limits of 3-7kpc provided by K05. 
Though suggestive of a disk-borne cluster, we find a non- 
negligible probability (10%) that a globular cluster passing 

through the disk would have this velocity. 

Our analysis of the Red Clump detected bv llvanov et al.l 

(2005) shows again that, within the accepted distance range, 
the cluster could have an age between several x lOOMyr 
and lOGyr. However, ages older than lOGyr appear to be 
ruled out as this would place the cluster in front of the 
nearby molecular cloud, incompatible with the cluster's ex- 
tinction. This is one argument against the cluster being an 

2 In calculating we have assumed that the ratio of mass to V - 
band luminosity M/Ly ~ 2 {McLaughlin fc van der Marej2 005l. 
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Figure 7. The if -band mass-to-light (M/L) ratio of a cluster 
as a function of age for three different m etalli cities, using the 
calcul ations of iMaraston fe Thoma and IMaraston et al.l 

j2003h . As the M/L we measure for GC01 depends on the cluster's 
physical size, which in turn depends on the distance, we indicate 
the cluster's full possible M/L range with the shaded area. The 
maximum M/L when using K05's measure of the cluster size is 
indicated by the long-dashed line. The arrow indicates the lower 
limit to the cluster age from the Red Clump analysis in Sect. [3751 



old globular. In addition, the cluster's high density, as well 
as marginal evidence for rotation, are both suggestive of an 
age much younger that that typical of globular clusters. 

We attempted to measure the RGB slope for a sample of 
synthesized clusters with a range of a ges and meta l licities , 
in order to compare to the results of llvanov et al.l (2005). 
However, large uncertainties in the measured slope caused 
by the Red Clump stars prevented us from reaching any firm 
conclusions, since the observations were consistent with both 
scenarios for the cluster's age and metallicity. 

In order to examine the nature of GC01 further, in 
Fig. [7] we compare the cluster's mass-to-light ratio (M/L) 
to predictions made by the stell ar evolutionary model s 
of IMaraston fc Thomas! (|2000h and IMaraston et al.l (|2003h . 
The figure shows how M/L varies in a cluster with a Kroupa 
IMF as a function of age at different metallicities. Our mass 
measurement of GC01 depends upon the cluster's size, and 
so depends linearly on the cluster's distance d; while the 
cluster's luminosity goes as d 2 . Therefore, M/L oc 1/d. For 
L, we use the integrated 7-C-band luminosity of K05, which 
we checked by re-measuring from both 2MASS and UKIDSS 
images, being careful to treat properly the effects of satu- 
rated stars. 

In Fig. [7] we indicate the observed range of M/L for 
the upper and lower distance limits. We see that, for ages 
greater than ~2Gyr, there begins to be a large disagreement 
between the observed M/L and the model predictions. Even 
using the larger cluster size as measured by K05 (indicated 
by the dashed line in Fig. [7} the maximum allowed age is 
3Gyr. This is consistent with our analysis of the Red Clump, 
which suggested an age between 300Myr - lOGyr (see Sect. 
13 . 3 fl - Therefore, we find a cluster age of 0.3-2 Gyr, which 
strongly indicates that GC01 is not an old globular clus- 



ter. This agrees with our findings of extremely high stellar 
density and marginal evidence for rotation. 

If, rather than being an old globular, the cluster is in- 
stead part of the disk population, we can use take the kine- 
matic distance of 5.0±0.9kpc and use Fig. [5]to further con- 
strain the age of the cluster to between 0.4-0.8 Gyr, which is 
perfectly consistent with the M/L ratio. We can then also 
constrain the mass of the cluster to (8±2)xl0 4 M Q . 

GC01 therefore represents the first detection of a very 
massive (M > 1O 4 M ) intermediate-age (10 8 - 10 9 yr) clus- 
ter in our Galaxy. Since s tar clusters are gra dually dis- 
rupted over their lifetimes ijLamers et al.ll2005h it is likely 
that the initial mass of GC01 was much greater that indi- 
cated by its current velocity dispersion, possibly by an order 
of magnitude. This would make it by far the most massive 
cluster known to have formed in the last lOGyr. Interest- 
ingly, the age of the cluster roughly coincides with peaks in 
the star-formation rate of the Solar neighbourhood and of 
both Magellanic Clouds (MCs) approx i mate l y 400 Myr ago 
l|Lamers et al.ll2005l ; lHarris fc Zaritskvl [200j . 120091 ). It has 
been speculated that these starburst events were caused by 
an interaction between the Galaxy and the MCs, and it is 
possible that the formation of GC01 was triggered during 
this episode. 



4.1 Future studies 

There are some caveats to the M/L analysis. Firstly, though 
the cluster's M/L is anomolously low for a globular, if it had 
passed through the Galactic disk numerous times, and was 
mass segregated, one may expect the lower mass stars at 
large distances from the cluster core to have been stripped 
away, which would drive M/L downward. Also, M/L may 
be somewhat sensitive to the precise nature of the IMF at 
sub-Sol ar masses (for a rev iew of the IMF in globular clus- 
ters, see lBastian et al]|2010l ). and the if-band luminosities of 
RGB and AGB stars. For these reasons, the nature of GC01 
still open for debate. Below, we suggest further experimen- 
tal tests that should provide incontrovertible evidence as to 
the nature of GC01. 

Metallicity studies: A cluster with an age typical of 
globular clusters would be expected to have a metallic- 
ity well below Solar (lo gjZ/Zp,) between -1.5 and -0.5, see 
iBrodie fc Stradej2006l ). On the other hand, an intermediate 
age cluster's metallicity would be much closer to Solar. Anal- 
ysis of the cluster members' spectra should readily indicate 
which regime the cluster belongs to, since the strengths of 
metallic li nes in the infrared a re very sensitive to abundance 
levels (e.g. lDavies et al.ll2010T ). 

Deep photometry: In Fig. [4] we show that, with cur- 
rent photometric data, we are unable to distinguish be- 
tween the CMDs of the two regimes. However, deeper pho- 
tometry (K > 19) would detect the location of the main- 



3 We note that, while llvanov et"aT, I ll2005l l attempted to measure 
the metallicity of GC01 from the slope of the red giant brach, this 
technique is calibrated only for old globular clusters. Attempts are 
currently underw ay to generaliz e this relation for clusters with a 
range of ages, see Sh arma et al. I ll2010l') . 
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sequence turn-off (MSTO) and it's position relative to the 
Red Clump. Going deeper still (K > 20), one is able to de- 
tect the kink in the main-sequence caused by the onset of 
molecular hydrogen absorption, which in conjunction with 
the MSTO has been shown t o be a powerful d iagnostic of 
age, distance and metallicity (Bono et al„ 21111);. An analy- 
sis involving two-colour photometry would be preferential to 
simple analysis of the cluster's luminosity function, as dif- 
ferential extinction and field star contamination would make 
the luminosity function difficult to interpret. From compar- 
isons to evolutionary predictions such as those in Fig. [4] 
deep IR photometry should provide the best estimate of the 
cluster's age and distance. 



5 CONCLUSIONS 

Using high-resolution near-IR spectroscopy of over 50 stars 
in GLIMPSE-C01, we have derived a dynamical mass for 
the cluster of (8 ± 3)xlO 4 M . Using our observations in 
conjunction with those in the literature, we have attempted 
to determine the nature of the cluster: whether it is an old 
globular cluster passing through the disk, or an intermediate 
age disk-borne cluster. The cluster's radial velocity, which 
could have ruled out the intermediate age possibility, in- 
stead indicates that the cluster is co-moving with the disk. 
Our analysis of the cluster's Red Clump, the mass-to-light 
ratio, as well as the marginal evidence we find of cluster 
rotation, all indicate an age of < lGyr. From our results 
we conclude that the cluster is part of the disk population, 
and we use the kinematic distance to constrain the cluster's 
age to 400-800Myr. The cluster is therefore the Galaxy's 
most massive intermediate age cluster discovered to date. 
In addition to our observations, we have outlined future ob- 
servational investigations capable of providing a critical test 
our hypothesis. 
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